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ABSTRACT 



An error probability analysis is performed for an M - ary orthogonal frequency 
shift keying (MFSK) receiver employing fast frequency-hopped (FH) spread spectrum 
waveforms transmitted over a fading channel with partial-band interference. The 
partial-band interference is modeled as a Gaussian process. Wideband thermal noise is 
also included in the analysis. Diversity is performed using multiple hops per data bit. 
Each diversity reception is assumed to fade independently according to a nonselective 
Rician process. A nonlinear combination procedure referred to as self-normalization 
combining is employed by the receiver to minimize partial-band interference effects. 
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I. INTRODUCTION 



Spread spectrum communication systems are designed to use signals having a 
much wider bandwidth than ordinary communication systems. Some of advantages 
of spread spectrum signals are the energy gain achivable against a narrow band jam- 
mer, low probability of detection (LPD), low probability of intercept (LPI), multiple 
access operation, and the capability to overcome the effects of multipath [Ref. 1]. 
The primary disadvantages of spread spectrum signals are the time and frequency 
synchronization required by the receiver which make it difficult to implement the 
system and the wide frequency band requirement. There are three primary spread 
spectrum communication methods: time hopping (TH), frequency hopping (FH) and 
direct spreading (DS). This thesis presents an error probability analysis of two fast 
frequency-hopped M-ary orthogonal frequency-shift keyed (FFH/MFSK) systems, 
where the term ‘fast’ implies that the frequency-hopping rate is greater than the 
symbol rate. The first to be examined is a conventional FFH/MFSK system and, 
the second is a self-normalization FFH/MFSK system. Both systems use noncoher- 
ent detection. The channel is assumed to experience both fading and partial-band 
interference. The performance of the two receivers is compared with each other as 
well as with another type of FFH/MFSK receiver referred to as a noise-normalization 
receiver [Ref. 2]. 

A. Conventional Receiver and Self-Normalization Receiver 

In the conventional FFH/MFSK receiver, illustrated in Fig. 1.1, the input 
signal passes through a wideband radio frequency (RF) filter and is mixed with the 
frequencies which are created pseudo-randomly in the frequency synthesizer to dehop 
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Figure 1.1: Block Diagram of Conventional FFH/MFSK noncoherent Re- 
ceiver 
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Conventional FFH/MFSK Receiver 



and convert the signal frequency down to intermediate frequency (IF). The bandwidth 
of the RF filter (IT) must be wide enough to pass all possible hopping frequencies, 
while the bandwidth of the IF filter ( B ) needs pass only one hopping frequency. 
Thus, the bandwidth of the RF filter (IT) is approximately equal to N B, where N is 
number of the hopping frequencies. The receiver has M different branches for each of 
the M different symbols. The IF signal is demodulated with a bandpass filter and a 
quadratic detector. The output signal from each quadratic detector is sampled L times 
and combined to provide a decision statistic between the M available outputs. The 
transmitter uses fast frequency-hopping to provide L independent samples through 
frequency diversity. The receiver is assumed to dehop the signal with perfect timing. 

In the self-normalization receiver, illustrated in Fig. 1.2, each of the M quadratic 
detector outputs are divided by the sum of the quadratic detector outputs from all 
M channels. This normalization factor is the only difference from the conventional 
receiver. When partial-band interference present, the hops that are corrupted by the 
partial-band noise jamming dominate the performance of the conventional detector 
since each hop is equally weighted in the overall decision statistic. The normalization 
factor minimizes this effect in the self-normalization receiver. As with the conven- 
tional FFH/MFSK receiver, the L samples at the output of the M quadratic detectors 
are combined to provide a decision statistic between the M available quadratic detec- 
tor outputs, and the receiver is assumed to dehop the signal with perfect timing. 

B. Fading 

For both receivers investigated in this thesis, the received signal is assumed 
to be a combination of a single, nonfaded (direct) component and many reflected 
(diffuse) signals. By modelling the fading as Rician, in which the signal is considered 
to consist of a direct signal component and a diffuse signal component, we obtain a 
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Figure 1.2: Block Diagram of Self-Normalization FFH/MFSK noncoherent 
Receiver 
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Self-Normalization FFH/MFSK Receiver 



general result that is valid in the limit of large direct-to-diffuse signal power ratios 
for Gaussian channels and in the limit of small direct-to-diffuse signal power ratios 
for Rayleigh fading channels as well as the general case where the effects of both the 
direct and diffuse components of the signal must be included in the analysis. For a 
Rician faded signal with amplitude the probability density function is [Ref. 3] 

/.,(«) (ID 

where a 2 is the power of the direct signal component, 2cr 2 is the power of the diffuse 
signal component, */(•) is the unit step function, and /o(*) is the zeroth-order modified 
Bessel function of the first kind. When the direct signal component power is zero, 
which occurs in times of deep channel fading on the channel such as when the direct 
communication path is blocked by terrain or other obstacles, the probability density 
function reverts to the Rayleigh probability density function [Ref. 3] 

/*»(<**) = ^exp (~^) u ( a k) (1-2) 

Also in this thesis, we assume that the smallest spacing between frequency hop 
slots is larger than the coherence bandwidth of the channel [Ref. 4]. As a result, each 
hop of a symbol experiences independent fading. In addition, the channel for each 
hop of a symbol is modeled as a frequency-nonselective, slowly fading Rician process. 
This implies that the signal bandwidth is much smaller than the coherence bandwidth 
of the channel and that the hop duration is much smaller than the coherence time of 
the channel [Ref. 5]. The latter assumption is equivalent to requiring the hop rate 
to be large compared to the Doppler spread of the channel. As a result, the signal 
amplitude can be modeled as a Rician random variable that remains fixed at least 
for the duration of a single hop. 
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C. Partial-Band Interference 



One of the applications of frequency-hopped spread spectrum signals is to reduce 
receiver performance degradation due to narrowband interference (either intentional 
or otherwise). This is accomplished by transmitting each symbol at a different carrier 
frequency according to some apparently random pattern known only by the trans- 
mitter and the receiver. The bandwidth of the frequency-hopped signal is generally 
much larger than the signal bandwidth in the absence of frequency-hopping. Conse- 
quently, in order to interfere with the entire spread spectrum band, given a fixed total 
interference power the noise power spectral density of the narrowband interference is 
reduced. Rather than interfere with the entire spread-spectrum bandwidth simultane- 
ously with the consequent reduction in noise power spectral density, the narrowband 
interference may affect only a portion of the total spread spectrum bandwidth at 
any one time, and the portion of the total spread spectrum bandwidth experiencing 
narrowband interference may also change in an apparently random manner. This 
is referred to as partial-band interference. In the case of partial-band interference, 
only some of the transmitted symbols will be affected by narrowband interference, 
and the question arises as to whether performance can be improved by implementing 
diversity in the form of fast frequency-hopping; that is, each information symbol is 
transmitted multiple times, and each transmission is at a different apparently random 
carrier frequency within the overall spread spectrum bandwidth. 

As mentioned previously, in this thesis the performance of two different fast 
frequency-hopped M-ary orthogonal frequency-shift keyed (FFH/MFSK) noncoher- 
ent receivers is investigated. The performance of the conventional FFH/MFSK non- 
coherent receiver in a Rician fading channel with only wideband noise (no partial- 
band interference) has been previously investigated [Ref. 6]. For the self-normalized 
FFH/MFSK receiver, system performance is obtained as an upper bound on the 
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probability of bit error, while for the conventional FFII/MFSK receiver, the actual 
probability of bit error is obtained. For both systems, the transmitter is assumed to 
implement 

In addition to wideband Gaussian noise, the channel is assumed to be affected 
by narrowband noise in the form of partial-band interference. The interference that 
is considered in this thesis may be due to a partial-band jammer as well as other 
unintended narrowband interferences. The partial-band interference is modeled as 
additive Gaussian noise and is assumed to corrupt only a fraction 7, where 1 > 7 > 0 , 
of the entire spread spectrum bandwidth at any one time. This is illustrated in Fig. 
1.3 where 5 represents the signal power level, A r j/2 the jamming noise plus thermal 
noise power spectral density, A r o/2 the thermal noise power spectral density, B the 
bandwidth of one frequency cell, and W the entire frequency band. In addition, 
partial-band interference, when present, is assumed to be present in each branch of the 
MFSK demodulator, and the fraction of the spread spectrum bandwidth experiencing 
partial-band interference is assumed to be the same for all hops of a symbol. Hence, 
7 is the probability that partial-band interference is present in all M channels of the 
receiver, and 1 — 7 is the probability that partial-band interference is not present in 
all M channels of the receiver. If Nj /2 is defined to be the average power spectral 
density of narrowband interference over the entire spread spectrum bandwidth, then 
7~ 1 Aj/2 is the power spectral density of partial-band interference when it is present. 
Thermal noise and other wideband noise that corrupt the channel are modeled as 
additive white Gaussian noise, and the power spectral density of this wideband noise is 
defined as A r 0 /2. Thus, the power spectral density of the total noise is 7" 1 Aj/2 + A 7 0 /2 
when partial-band interference is present and Aq /2 otherwise. If the equivalent noise 
bandwidth of each bandpass filter in both the conventional and the self-normalized 
FFH/MFSK demodulators is B Hz, then for each hop the signal is received with noise 
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Figure 1.3: Thermal Noise and Partial-Band Jamming Noise Model 
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Thermal Noise and Partial-Band Jamming Noise Model ( 



of power N 0 B with probability 1— 7 when partial-band interference is not present and 
with noise of power (7 -1 Nj + N 0 )B with probability 7 when partial-band interference 
is present. 

The symbol rate is R s = R^j log 2 M where M is the order of the MFSK modula- 
tion and Rb is the bit rate. Since the MFSK signal is assumed to perform L hops per 
symbol, the hop rate is Rh = LR S . The equivalent noise bandwidth of the bandpass 
filters in each of the M channels of both the self-normalized FFH/MFSK receiver and 
the conventional FFH/MFSK receiver must be at least as wide as the hop rate, and 
in this thesis we use B — Rh- The overall system bandwidth is assumed to be very 
large compared to the hop rate. 

D. Fast Frequency-Hopping 

There are two types of frequency-hopping: fast frequency-hopping and slow- 
frequency hopping. In fast-frequency hopping the carrier frequency changes (hops) 
more than once per data symbol and all hops within a symbol duration are combined 
in the receiver to provide the sample statistic, but in slow frequency-hopping the 
carrier frequency changes only once per data symbol. Hence, the terms of ‘slow* 
and ‘fast* do not mean absolutely how fast the carrier frequency changes. We define 
diversity of order L as the case where each M-ary symbol is transmitted L times over 
the channel. Thus, if Eh is the energy for each trasmission then each M-ary symbol 
requires a total energy of 

E s = LEh ( 1 . 3 ) 

where E s is the symbol energy. The diversity factor is a key point of fast frequency- 
hopping communication systems, and the object of this thesis is to determine under 
what conditons fast frequency-hopping systems can provide immunity to partial-band 
interference. 
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II. ANALYSIS 



A. Determination of the Probability of Bit Error 
1. Conventional Receiver 

An analysis of the conventional FFH/MFSK receiver requires, to begin 
with, a description of the samples at the output of the M quadratic detectors be- 
fore diversity combining. The partial-band interference, when present, is assumed 
to be present in all channels. Assuming equally likely M - ary symbols, we have the 
probability of symbol error given L independent hops as 

= )V(i-7) t “'>,(i) (2.i) 

where P s (?) is the conditional symbol error probability given that i out of L hops 
are jammed. Assuming the signal is present in branch one, we have the probability 
density function of the random variable X\ k representing the output of the quadratic 
detector [Ref. 3]. 

fx lk {xik\ak) = ^2 ex P (- Tl *9^2 2a * ) 7 ° ( 2 . 2 ) 

where y/2 cik is the signal amplitude, a\ is noise power in hop fc, and /o( # ) is the 
zeroth-order modified Bessel function of the first kind. The unconditional probability 
density function is found by evaluating 



fx lk {xik) = / fx lk (xik\ 0 'k)fA k (ak)dak (2.3) 

jo 

The amplitude of the signal is modeled as a Rician distributed random variable with 
a probability density function as given by (LI). Substituting (1.1) and (2.2) into 
(2.3), we have 



fx lk {xik) = 



1 

2 K 2 + 2a2) 



exp 



%\k + 2a^. 

'2(^ + 2^). 



Io 



cik \/2xj^ 

.ol + 2a* 



u(x lk ) 



(2.4) 
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The probability density function of each of the M — 1 channels with only Gaussian 



noise is 



f Xmk ) = 2~2 eX P 



%mk 

"M, 



I u(Xyk) 



(2.5) 



which is easily obtained by replacing the signal amplitude a* with 0 and replacing Xu 
with where m = 2, 3, 4... .M . Next, the L samples are combined to obtain the 

final sample statistic. Assuming i out of L hops are jammed, we obtain the decision 
variable for branch 1 



A'i 



EAu- 

fc=i 
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(2.6) 



where the superscript ^ implies a jammed hop and ^ implies a non-jammed hop. 
Since each hop is assumed to be independent, we have 



fx,( x i) 



( 8 ) 



f x m(x S?) 

Ik 






where is the probability density function when jamming is present, f y (2)(x\V) 

when no jamming is present, ® implies convolution, <g>? implies an i — fold convolu- 
tion, and ®(L — i) implies an ( L — i)-fold convolution. Similarly, for the noise only 
branches, we have the random variables 
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and the probability density function 
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For MFSK. the probability of symbol error is [Ref. 5] 



pm = i - /' M*i\‘)\r 
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Equation (2.10) requires fx 1 (xi|i) and fx m { x m )• In order to obtain fx 1 (xi |i), 



we first rewrite (2.7) 

/jri(*iio = /x 1 (*ii*)®yx 1 (*ii*) 

Two probability density functions are computed separately 



( 2 . 11 ) 
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where F\ ik (s) is Laplace transform of fx lk (xik\i) and C 1 is the inverse Laplace 
transform operation. 
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(2.13) 



1 + 2s(<j\j + 2<j 2 )_ 

where crjL represents jamming noise power plus thermal noise power. The inverse 
Laplace transform of (2.13) yields 
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«• + 2^ 2 )J 


I t-l 





(2.14) 



Similarly 



fxA x il*) = ex P 



2 (L — i)a 2 + Xi 


'x t (al + 2a 2 )' 


1 

(M 

b 

CM 

+ 

IN-* 

^b 

i 


• <s> 

1 



X 



1 


L-t+ 1 
2 

Tr . 


\J2{L - i)a 2 k x x 


2(a\ + 2a 2 ). 


-* L—i — 1 


a\ + 2a 2 



(2.15) 



where <j\ represents thermal noise power. In order to obtain fx 1 (x\\i ) 1 (2.14) and 
(2.15) must be convolved. In general, this cannot be done analytically. 
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For the noise branch, the probability density function is 



2< 



fXm{ X m) — 

The Laplace transform of (2.1G) is 

FxM = 

where 



1 2a l 

e k > 



0i 



(8) 



2L e -s* 

2 °l 



L—i 



®(L- 1) 



(2.16) 



d‘> 1 


l 


l 

1 


1 

+ 

l 




_s + 4 2) . 



/9 (1) - — /? (2) - — 

Pk ~ 2 < ’ Pk ^ l 



(2.17) 



(2.18) 



Equation (2.17) can be expanded in a Heaviside partial fraction expansion to obtain 

K u 



Fx m (s) = 



+ 



An A 12 

* + d" ( S +tf>) 2 



+ ••• + 



A 2i A 



22 



5 + 4 2) ' (5 + 4 2) ) 2 



+ ••• + 



(s + tf ) )’ 

•^2(L-i) 

+ d 21 )' 1 '"’ 



(2.19) 



where the coefficients K\j, j — 1 , 2 , •••,* and A'n- , A' = 1 , 2 , — i depend on 

A and i. Specific coefficients for several values of L and i are given in Table 2.1. 
Applying the inverse Laplace transform to (2.19), we obtain 



.fx m (Xm) — £ Xm ( All + A 12 X m + 






+e~ p k 



P[ 7) Xm 1 L' 



A 21 + A 22 X m + 



2 ! 

A'23* 



+ • • • + A 1,-7 



X 



1-1 



(»*-!)!, 



2 ! 



+ A 2(L_t)‘ 



„L-t-l 



(2.20) 



(A-/-l)!y 

We can now compute the second integral in (2.10) using (2.20) with the 



result 
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L, i 


Coefficients 


II 

1 — 1 
os . 

II 

o 


A 2 1 = W 


L = 1, i = 1 


An = W 


L — 2, i = 0 


A'j, = 0 

= (4 21 ) 2 


L = 2, i = 1 


qK 1 ) ft( 2 ) 

K u = -4/* 12) 

Zv 21 = — Z\ii 


II 

to 

os . 

II 

to 


A'„ = 0 

a-,2 = (tf ] y 


II 

5° 

<*o . 

II 

o 


I<n = 0 

I<22 = o 

I< 23 = (4 2) ) 3 


Z, = 3, i = 1 


Ki — 3^ ( ^ 2> ^ 


Pk v-^>+/?‘ 2 v 

A ” “ 

Am = /W - A’n/Sf 1 - fS 


II 

OS . 

II 

to 


Zv> — ( ^ Y R^ 


V-z^+^V fik 

12 “ -/Kef 

A'l. = PlM'-KiA" - fg 


II 

co 

OS . 

II 

CO 


Zvn = 0 
Ku = 0 

A'i3 = (4 1J ) 3 



TABLE 2.1: Coefficients of Heaviside Partial Fraction Expansion 
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L - £ { p<» y 



- V A, T — 

h " sf-i-'iwr 1 



P= 1 



i — 1 — r 



+ E 



L-i R 



2 p 



s (tfV 



P 
L-i 



5 A ^ e ’ S (L - i - 1 - r)!(4 2, )'+> 

So the final form of P 5 (z) which must be evaluated numerically is 

2ia\ + a:i \ 



p= i 



L — i — l—r 



p - (i) = >-r[M- 



2( CT^- +2tr 2 ) / 



X 



*i« + 2<r 2 )' 



1M 

tal 



i-l 

2 



1 



1±1 

2 



2(<r? +2<r 3 )J 



exp - 



fc / V 

2(Z, - i)a| + a-! \ ( (erf + 2cr 2 ) 






2(al + 2a 2 ) 



(L - i)al 



X 



X 



,2W + 20, 



, /" \/2(A- | 

t — ' ( <r? + 2 <t 3 J | 



P =1 



E^t-E E r. 



i — l—r 



P= 1 (/^Jt ) P P=1 



'„(i-l -.')!(/?<■ V 



L-i r- L~t , 0 v 

+ E - - E av* '■ 



(/5f , ) p „.i 

P=1 



X 



E 



_L— i— 1— r 



o (L — « — 1 — r)!(/?{ 2 >) r+1 



A/-1 






(2.21) 



(2.22) 



2. Self-normalization Receiver 

The difference between the self-normalization receiver and the conventional 
receiver is the self-normalization factor which is the sum of all M channel outputs 
after sampling at the outlet of the quadratic detectors. The new decision statistic is 
the combination of self-normalized quadratic detector outputs. In order to find the 
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probability density function of the random variables representing the self-normalized 
quadratic detector outputs, we first define the random variable 



Vk — (X 2k + Xzk + • * • + X\4k) ” Xmk , (m 7 ^ 1) 



(2.23) 



The probability density function of the sum of independent random variables is the 
convolution of the probability density functions of those random variables. Thus, the 
probability density function of the random variable V k requires M — 2 convolutions: 



fv k (v k) = fx 2k (x 2 k) 0 fx ak (x 3k ) ® ® fx, k (Xik) ® 

® fx Mk (x M k) {i + rn) 



(2.24) 



Assuming that all the probability density functions of the branches containing only 



noise are the same, we have 



fv k (v k ) = [fx, k (xik)]® {M 2) (i + 1 or m) 



(2.25) 



The Laplace transform of fx, k {xik ) is 



Hence 



= 



Fv„(s) = 



2 <7?S + 1 



1 



M —2 



2(j?5 T 1 



Taking the inverse Laplace transform of (2.27), we get 



fv k (vk) 



M ~ 2 \M - 3 



= Mr M 



exp 









(2.26) 



(2.27) 



(2.28) 



M) 

As can be seen from Fig. 1.2, the random variables that represent the self- 
normalized outputs of each of the M channels after diversity combining are obtained 
as 

L 

C4 = £tW = l,2,...,M (2.29) 

jt=i 
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where U{ k , i = 1,2, . . . , A/, A’ = 1,2, ...,L are the random variables that repre- 
sent the self-normalized outputs of each of the M channels for hop k of a symbol. 
These random variables are expressed in terms of the independent random variables 
Xik , i = 1,2,..., M, k = 1,2,...,// representing the M independent outputs of the 
M quadratic detectors for each channel of the receiver for hop k of a symbol as 



u ik = 



x tk 



(2.30) 



Xik + X 2k + . . . + X\ik 
where 0 < U lk < 1 since 0 < X xk < oo ,i = 1,2, Using 2.29, we can express 

Pr(Ui < U m \i) as 

Pr {U, < U m \i ) = Pr (jr u lk < jr = Pr (£ Un - U mk < 0|A (2.31) 

\A-=I k=l / \k= I / 

Defining the random variable = Uik — U m k , we obtain 



Pr {U\ < U m \i) = Pr^ Zi*<0|*') 



= Pr [Z\ < 0|f) 



(2.32) 



where 



L 



z^ = Y, Zik 

k = i 



is our alternative decision variable. 

From 2.30 and the definition of Zu- above, we have 



(2.33) 



Zk = U\k ~ U m k 

Au- ~ Xmk 

Au- + X2 k + • • • + A An- 
A ik — Xmk 
X u- + X mk + V k 



(2.34) 



where —1 < Zk < 1. This random variable represents the difference between the 
power in the branch containing the signal and the power of one noise only branch 
which are detected by quadratic detectors and normalized by the self-normalization 
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factor. We use the auxiliary variable method [Ref. 7] to obtain the probability density 
function of the new random variable Zi k , where we define the two auxiliary variables 



W k = 14, Y k = X lk 



(2.35) 



Reforming (2.34) by replacing V k and X\ k with (2.35), we have 

Y k { 1 - Z lk ) - Z lk W k 



X mk — 



1 + Z\ k 



(2.36) 



and the Jacobian of the transformation is 



(1 +*u) 2 

2 yk + w k 



(2.37) 



Now the joint probability density function for the random variables Z ik , W k , and 
Y k is obtained in terms of the joint probability density function of the independent 
random variables X\ k , X mk ■ and V k as 



fz lk Y k w k {z lk ,y k ,w k \i) 



jj|/A' lfc (yr|0/x mfc 



'yk( 1 - zik) ~ z lk w k \ 

, TT^ j Mwi) 



tyk + W k 

(i + Zlk y 



X 

X 

X 



1 / Vk + 2 a\ \ 

2{<J 2 k + 2cr 2 ) eXP V 2(cr^ + 2cr 2 ) ) 
1 { yfc(l ~ z lk ) + 

2 o\ 6XP V 2a^(l + zn) J 
( 1 \ M ~ 2 (w k ) M ~ 3 ( w k 

\2al) (M — 3)! GXP \ 2a 2 k 



I o 



/ aW2yk_ \ 

W + 2* 2 ) 



(2.38) 



Before evaluating the probability density function of the random variable Z \ k , we must 
be careful in handling the range of the variables. Note that the range of the variable 
X m k is — co < X m k < oo, but recall that fx mk (xmk) is nonzero only for the range 
Xmk > 0. For the range of — 1 < Z\ k < 0, X m k is always greater than zero; but 
for the range of 0 < Z\ k < 1, the range of X m k is — oo < X mk < oo. Consequently, 



18 



special care must be taken in evaluating fz lk ( z ifc)- The probability density function 
of the random variable Z\ k is obtained from 

n OO | 

jj\fz lk Y k w k (zn,yk,Wk\i)dw k dyk (2.39) 

when — 1 < Z\k < 0- Since fx mk { •) is nonzero only for positive values of its argument, 
then 



n ja-(i-zu)/*ifc 1 

Zi k Y k w k ( z ik, Vk, Wk)dwkdy k (2.40) 

when 0 < Zu- < 1. Evaluating the integrals in (2.39) and (2.40), we have the 
probability density function of random variable Z\ k conditional on i 



fz lk (zik\i) = 



M-2 



r / -p fc (i - zjk) A 

2 + ^(1 — z\k) ) ^ \2 + £fc(l — z\k)) 



(1 + Z\k) 



(m-9 + 2(1+ . 6) .. 

\ + 2 + 6(1 ~zn) 

-g{zik)u{2 lk ) 



l + 



Pk{ 1 + Zik) 



(i + 6(i-^u))(i + 6) 



(2.41) 



where 



g(zn) = 



'U-( 1 — Z\k) 1 2 (1 + ik) 



\M-2 



(1 + 6.(1 — Z\k )) fl/ ! (1 + Zik) ^ + 6(1 — z \k ) , 

( i / \ M-2 

I 1 ( PkZ\k \ 

X 



— pk{ 1 — Zik ) 



1 / pk-z\k y 

(M- 3)! \,(1-6(1-*i*))(1+6)J 



M —3 ( 

+ E 
m = 0 t 



+ 



1 + Z\k 



(M - 2)(M - 2)! 

(M-2 - m)!(77i!) 2 

pr-~r- 

, 0 + 6(1 — ^ifc))(l + 6) 

2Z\k 

{1 - z lk )(M - Z - m)\ [(1 +6(1 -*i*))(l + 6) 



, 1 + 6( 1 — Zik) t 

PkZlk 



M — 2 — m 



M —3—m 



+ E 



71/ — 3 — m 



PkZ\k 



X 



- (M -3-m- ,)!(,!)= VO + <i(l - 2ii))(l + fO/ 
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X 



2z lk (M -2 -m) 2 
(1 - z xk ){M —2 — m — q) 

(1 + ~ g ifc))(l + Z\k){M — 2) 

(1 - «ifc)(l + ik) 



and 



(2.42) 



Pk 

ik 




(2.43) 



are the signal-to-noise ratios of the direct signal component and the diffuse (or faded) 
signal component, respectively. 

Having obtained the probability density function of the random variable 
Ziki we need to consider diversity to get the overall output sample statistics. The 
random variable Z\ is defined 



Zi = Ezi k 



k=l 



= E zJl> + L zl? 

k = t + 1 



(2.44) 



k—i 



and the probability density function of the random variable Z\ is obtained as 



- -i Qi r 

fzi (*i 10 = f z o)( z n) <g> [fz™ (4* ) 



-i ®(L-i) 



(2.45) 



It is difficult to evaluate equation (2.45), so we leave this problem for numerical 
analysis. 

Now the probability of symbol error of the self- normalized FFH/MFSK 

receiver is 



P s (i) = Pr(U \ < U 2 (J < tO U •••U < U mV) (2-46) 

For M > 2, an analytic solution of (2.46) is not possible because the random variables 
Uii ^3? * * * Ufrf—i are not independent. Hence, we use the union bound to obtain 
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an upper bound on system performance. From (2.46) 

Ps(>) < Pr(U r < U 2 \i) + Pr(U, < U 3 \i) + ••■ + Pr(Ux < U M \i) (2.47) 

Since we assume the sample statistics of all noise only branches are the same, then 
(2.47) becomes 

Ps(i) < (M - 1) Pr(U x < U m \i) (2.48) 

where m — 2, 3, 4, • • • , M . Substituting (2.32) into (2.48), we have 



P s (i) < (M - 1 )Pr(Zx < 0|t) 



(2.49) 



Since 



then from (2.49) 



Pr{Zx < 0) = f fzA z i \i)dz\ 

J — Lj 

Ps(i ) < (M — 1 ) /° fzA*\\')dzx 

J — L 



(2.50) 



(2.51) 



The probability of symbol error is related to the probability of bit error by [Ref. 5] 



9^-i 

P b (i) = P s (i ) x ,fc = log 2 M 



- P,(i) x 



2 ^ — 1 

M 



2 (M - 1) 

Substituting (2.48) into (2.52), we get 

M f° 

W) < — }_ L hM\i)i ■ 

Equation (2.53) must be evaluated numerically. 



(2.52) 



(2.53) 
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III. NUMERICAL ANALYSIS 



Results of the probability of error analysis for both the conventional receiver 
and the self-nomalization receiver are obtained for L = 1 to 6 hops per bit for several 
signal-to-noise ratios and direct-to-diffuse component ratios as a function of signal-to- 
interference ratio to provide an overall comparison of channel and interference effects 
on system performance for both receivers. In all cases the evaluation of either (2.22) 
or (2.41) is accomplished via a numerical integration routine. We consider L + 1 
different cases for a diversity of L hops. For example, for L = 1 the two cases to be 
considered are one non-jammed hop and one jammed hop. For L = 2 there are three 
cases: neither hop jammed, one non-jammed hop and one jammed hop, and both hops 
jammed. Additionally, the effect of the jamming ratio 7 on system performance is 
examined. We can relate the signal-to-noise ratio to the signal energy density-to-noise 
power ratio for a jammed hop as 

a\ + 2a 2 E b B log 2 M 



J kj 



No + NjhRs L 



(3.1) 



where 



a 2 : direct signal power component 

2er 2 : diffuse signal power component 

a l + 2< t 2 : total signal power 



u kj 

E b 

No 

n 3 

7 



thermal noise plus jamming noise power 
average signal bit energy density 
thermal noise power spectral density 
jamming noise power spectral density 
jamming ratio 
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B 



: IF bandwidth 



R s : symbol rate 

M : number of symbols in alphabet 

L : number of hops in a symbol 

For hops with no partial-band jamming 

a\ + 2<t 2 Eb B log 2 M 

W k ” W 0 W s L 

where cr 2 is the thermal noise power. Assuming the IF bandwidth ( B ) is equal to the 
symbol rate, we can rewrite (3.1) and (3.2) 

Eb log 2 M 



(3.2) 



a 2 k + 2 < t 2 



a h 



No + Nj/i L 



and 



(3.3) 



+ 2 <r 2 


E b log 2 M 
No L 


(3.4) 


For computational purposed, define 




= 1 


(3.5) 




Eb 

No 


= V 


(3.6) 




E b 

Nj 


= V 


(3.7) 




Q l 

2 a 2 


- *l> 


(3.8) 


log 2 M 


= cu 


(3.9) 


Substituting (3.5) through (3.9) into (3.3) and (3.4), we have 




2 a 2 = 




7 / 72 / u 


(3.10) 


(r; + 7 J /){%!> + 1) L 


rv 2 






(3.11) 


a k — 


(V + 7' / )(V’ + 1 ) L 






72 / 


(3.12) 


a k — 


(i) + 7 i/) 
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Substituting (3.10), (3.11), and (3.12) into the probability density function of the 
conventional receiver and the probability density function of self-normalization re- 
ceiver and using an fast Fourier transform (FFT) routine with a sequence length of 
N = 1024, numerical sequences of the probability density functions are generated for 
various signal- to-noise ratios of the jammed and non -jammed cases. 



A. Conventional Receiver 

Except for a single convolution, we have obtained analytic solution for the prob- 
ability density function of the decision statistic for the conventional FFH/MFSK de- 
tector. Unfortunately, the modified Bessel function terms make it difficult to obtain 
a complete analytic solution for the general case. Consequently, we use an FFT to 
perform the convolution numerically. Even though we have an analytic expression for 
the probability density function of the noise branch statistic (the last part of (2.22)), 
we use numerical methods to correspond with the signal statistic part which must be 
evaluated numerically. Hence, we use different equations for the numerical procedure. 
Equation (2.16) can be rewritten 

fx m (Xm) = fxl( x m) ® fxl( X m) (3.13) 



where 



II 

"7 

— $ 


1 

e *} 

[2 *ij 


II 


1 ,-%l] 

~ 7Z e k 

l 2a k J 









Using the Laplace transform, we get 

4 ; Iw = 

ffj>) = 



p, 



( 1 ) 1 * 



s + 






> + 



L-i 



(3.14) 

(3.15) 



(3.16) 

(3.17) 
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where 




1 



(3.18) 





1 



(3.19) 




The inverse Laplace transform of (3.16) and (3.17) yields 





(3.20) 




L—i 





m 

(L — i — 1)! 



(3.21) 



Now we have a different form of the probability density function of the random vari- 
able X m 



Instead of (2.21) we substitute (3.22) into (2.22) and analyze (2.22) numerically be- 
cause we need a numerical sequence of this formula to coincide with the numerical 
sequence obtained for (2.11) The length of the basic sequence is N=1024, and prior 
to convolution it is zero-padded to create a sequence of total length N=204S since the 
sequence is to be convolved one time. 

Numerical problems are encountered since as the value of the variable .Ti in- 
creases the exponential terms approach zero and the modifiied Bessel function terms 
approach infinity. Hence, the multiplication of the exponential term and the modified 
Bessel function term is invalid when the argument of the exponent is less than -745 
or the argument of the modified Bessel function is greater than 705 when 386Mat- 
lab is used. For example, in the situation described above, the value of exp(-745) is 
4.9470e-324, but the value of exp(-74G) is set to 0 by the computer. The value of 





L-i 



(3.22) 
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Bessel(705), where the 386 MATLAB function for the modified Bessel function 7o( •) 
is Bessel(#), is 2.2621e+304 and Bessel(706) is not defined. Even though the product 
of the two is either zero or an invalid number, the true value of a multiplication is a 
valid number. When this situation is encountered, the program evaluates the situa- 
tion and uses either zero or one (because the cumulative probability density function 
approaches one as the value of the variable increases). 

When two probability density functions are convolved by an FFT, if the valid 
range of one probability density function is much wider than the range of the other, the 
result is almost the same as the probability density function that has wider range. 
We must check for this situation whenever the parameters of the equation change 
while the program is running, skip the numerical convolution steps, and instead use 
the probability density function that has wider range as the convolution output. 

B. Self-Normalization Receiver 

The probability density function of self-normalization receiver (2.41) is relatively 
complicated and we must use an L-fold convolution to evaluate (2.44); an z-fold con- 
volution for jammed hops, an ( L — z)-fold convolution for nonjammed hops, and each 
result must be convolved. To implement the convolution numerically, the numerical 
sequences of the probability density function of the self-normalization receiver must 
be zero-padded to a total length of TV x L before the FFT is performed in order 
to preserve the desired linear convolution of (2.44). The transformed sequences are 
multiplied point by point to produce the desired seqence of 

FzM = K!wF X KlM] L " (3.23) 

where F is the Fourier transform of /. This result is then inverse fast Fourier trans- 
formed to yield the desired sequence for fz^zi). The procedures for L= 3 and L— 5 
are similar to that for L= 4 and L= 6, respectively, because the length of the sequence 
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for the FFT should be a power of two, and the sequences should be zero-padded to 
A r =4096 (2 12 ) for L = 3 and L = 4 and to A r =8192 (2 13 ) for L = 5 and L— 6. 

Equation (2.42) has many summations inside which must be carefully handled 
since numerical errors are likely when the value of the variable Z\k approaches one 
When z\k approaches unity, the term 1/(1 — z\k) approaches infinity. Therefore, we 
must group the terms properly to minimize numerical errors. Equation (2.42) is 
rewritten 



s(-u-) = 



•?H-(1 ~ ^u) A/ ~ 2 (l + jk) M ~ 2 
(1 + < ffc ( 1 — 2 U -)) M_1 (1 + z \k) 

(M - 2)(M - 2)! 

X m^o ( M - 2 - ™)\{m\y 

A/-3 



exp 



Pk ( 1 ~ Z\k) 

1 + 6( 1 ~ Zik) 

PkZ\k 



(i +6(1 - «1*))(1 +6) 



+ £ ex p 



m= 0 



~Pk{ 1 ~ Zlk) 

1 + 6(1 — Z\k) 



>?*((! - *u)(l +6)) A/_3 " m (l + 6)U + z lk ) m - 1 (M-2 -m) 



(1+^(1 -2 U -)) m_1 ' 



A / — 2 — m 



X 



£ 



M 



m 



<7 

M-3 



% (M - 2 - m - q)l(ql) 2 



PkZ\k 



(1 + 60 — -^U-))(l + ik) m 



+ £ exp 

m — () 

A/ — 3— m 



-Pk{ 1 - Zn) 
[1 + 6(1 -zik) 
i XI — 3 — m 



((1 - 2 l/c)0 + ik)) 



M — 3 — m 



z\k{ 1 + zi k ) 7n (M — 2) 



(1 +f(l -z\k)) M ~ 2 - 



5 (M - 3 - m - q)\(q\) 2 



PkZn 



.(1 + 6(1 — *ifc))(l + 6) 



(3.24) 



in order to minimize computational errors. 
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IV. RESULTS 



A FFH/MFSK receiver design that demonstrates strong immunity to the partial- 
band interference is the noise normalization receiver. This receiver requires a measure- 
ment of the noise power per hop, but it is not easy to measure the noise power precisely 
in practical situations. This measurement is not required for the self-normalization 
receiver; hence, it is easier to implement the self-normalization receiver. As will 
be seen, partial-band jamming strongly affects the performance of the conventional 
FFH/MFSK receiver, while the self-normalization FFH/MFSK receiver performs well 
against partial-band jamming. Probability of bit error analyses for both the conven- 
tional system and the self-normalization system assume worst case partial-band jam- 
ming. The channel is modeled for a moderately strong direct-to-diffuse signal ratio 
(a 2 /2a 2 = 10 ) which we will refer to as Rician, and as a primarily Rayleigh faded 
channel (a 2 /2a 2 = 0.01). For Rician fading the worst case jamming ratio typically 
increases as interference power increases. For Rayleigh fading worst case jamming 
corresponds to broadband jamming (jamming ratio 7 = 1 ) for all levels of jamming 
power. 

A. Performance of Conventional FFH/MFSK Receiver 

Graphs of the probability of bit error, including worst case, as a function of 
signal-to-jamming ratio for fixed signal-to-thermal noise ratios are obtained for 1 
through 6 hops per symbol. These results are obtained by numerically evaluating 
(2.44) for values of the jamming ratio ranging from 1.0 to 0.001 and retaining worst 
case performance for each value of Ei/N 3 . The results for a Rician fading channel 
are shown in Fig. 4.1 through Fig. 4.6. 
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Figure 4.1: The Performance of the Conventional FFH/MFSK nonco- 
herent Receiver at Signal-to-Thermal Noise Ratio = 13.35 dB, Hopping 
Number per Symbol = 1, Order of Modulation = 4 and Direct-to-Diffuse 
Signal Power Ratio = 10 
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Figure 4.2: The Performance of the Conventional FFH/MFSK nonco- 
herent Receiver at Signal-to-Thermal Noise Ratio = 13.35 dB, Hopping 
Number per Symbol = 2, Order of Modulation = 4 and Direct-to-Diffuse 
Signal Power Ratio = 10 
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Figure 4.3: The Performance of the Conventional FFH/MFSK nonco- 
herent Receiver at Signal-to-Thermal Noise Ratio = 13.35 dB, Hopping 
Number per Symbol = 3, Order of Modulation = 4 and Direct-to-Diffuse 
Signal Power Ratio = 10 
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Figure 4.4: The Performance of the Conventional FFH/MFSK nonco- 
herent Receiver at Signal-to-Thermal Noise Ratio = 13.35 dB, Hopping 
Number per Symbol = 4, Order of Modulation = 4 and Direct-to-Diffuse 
Signal Power Ratio = 10 
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Figure 4 . 5 : The Performance of the Conventional FFH/MFSK nonco- 
herent Receiver at Signal-to-Thermal Noise Ratio = 13.35 dB, Hopping 
Number per Symbol = 5 , Order of Modulation = 4 and Direct-to-Diffuse 
Signal Power Ratio = 10 
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Conventional FH/MFSK Receiver Performance 
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Figure 4.6: The Performance of the Conventional FFH/MFSK nonco- 
herent Receiver at Signal-to-Thermal Noise Ratio = 13.35 dB, Hopping 
Number per Symbol = 6, Order of Modulation = 4 and Direct-to-Diffuse 
Signal Power Ratio = 10 
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The signal-to-noise ratio, direct-to-diffuse ratio, and the modulation order are 
chosen to enable direct comparison with the results obtained in [Ref. 2 ]. We see that 
for a Rician channel that diversity is not effective in reducing the effects of worst case 
partial-band jamming. Worst case performance is poorer when diversity is used than 
for no diversity. As diversity increases, only the asymtotic value for high Et>/Nj ap- 
proaches the performance of the noise-nomalization receiver. The order of modulation 
has no effect with respect to providing immunity to partial-band interference. This 
is illustrated in Fig. 4.4 and Fig. 4.7 through Fig. 4.9. Generally both diversity and 
the order of modulatin are effective in improving system performance for full-band 
jamming (7 = 1 ) as shown Fig. 4.10 and Fig. 4.11. 

Generally both diversity and the order of modulation for a Rayleigh fading 
channel are effective in improveing system performance for full band jamming ( 7 = 1 ) 
as shown in Fig. 4.10 and Fig. 4.11. However, System performance, shown in Fig. 
4.12 through Fig. 4.16, is not improved by either diversity or modulation order when 
worst case partial-band jamming is assumed. When the signal-to-thermal noise 
ratio increases, the performance of the conventional receiver for worst case partial- 
band jamming is unchanged from that of lower signal-to-thermal noise cases except 
for the asymtotic limit of negligible E^/Nj. 
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Conventional FH/MFSK Receiver Performance 
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Figure 4.7: The Performance of the Conventional FFH/MFSK nonco- 
herent Receiver at Signal-to-Thermal Noise Ratio = 13.35 dB, Hopping 
Number per Symbol = 4, Order of Modulation = 8 and Direct-to-Diffuse 
Signal Power Ratio = 10 
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Figure 4.8: The Performance of the Conventional FFH/MFSK nonco- 
herent Receiver at Signal-to-Thermal Noise Ratio = 13.35 dB, Hopping 
Number per Symbol = 4, Order of Modulation = 16 and Direct-to-Diffuse 
Signal Power Ratio = 10 
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Figure 4.9: The Performance of the Conventional FFH/MFSK nonco- 
herent Receiver at Signal-to-Thermal Noise Ratio = 13.35 dB, Hopping 
Number per Symbol = 4, Order of Modulation = 32 and Direct-to-Diffuse 
Signal Power Ratio = 10 
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Conventional FH/MFSK Receiver Performance 




Figure 4.10: The Performance of the Conventional FFH/MFSK noncoher- 
ent Receiver at Signal-to-Thermal Noise Ratio = 20.0dB, Jamming Ratio 
= 1, Order of Modulation = 8 and Direct-to-Diffuse Signal Power Ratio 
= 0.01 
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Figure 4.11: The Performance of the Conventional FFH/MFSK noncoher- 
ent Receiver at Signal-to-Thermal Noise Ratio = 20.0dB, Hopping Num- 
ber per Symbol = 4, Jamming Ratio = 1 and Direct-to-Diffuse Signal 
Power Ratio = 0.01 
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Figure 4.12: The Performance of the Conventional FFH/MFSK nonco- 
herent Receiver at Signal-to-Thermal Noise Ratio = 13.35 dB, Hopping 
Number per Symbol = 1, Order of Modulation = 4 and Direct-to-Diffuse 
Signal Power Ratio = 0.01 
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Figure 4.13: The Performance of the Conventional FFH/MFSK nonco- 
herent Receiver at Signal-to-Thermal Noise Ratio = 13.35 dB, Hopping 
Number per Symbol = 3, Order of Modulation = 4 and Direct-to-Diffuse 
Signal Power Ratio = 0.01 
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Figure 4 . 14 : The Performance of the Conventional FFH/MFSK nonco- 
herent Receiver at Signal-to-Thermal Noise Ratio = 13.35 dB, Hopping 
Number per Symbol = 5 , Order of Modulation — 4 and Direct-to-Diffuse 
Signal Power Ratio = 0.01 
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Figure 4.15: The Performance of the Conventional FFH/MFSK nonco- 
herent Receiver at Signal-to-Thermal Noise Ratio = 13.35 dB, Hopping 
Number per Symbol = 4, Order of Modulation = 16 and Direct-to-Diffuse 
Signal Power Ratio = 0.01 
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Figure 4 . 16 : The Performance of the Conventional FFH/MFSK nonco- 
herent Receiver at Signal-to-Thermal Noise Ratio = 13.35 dB, Hopping 
Number per Symbol = 4 , Order of Modulation = 32 and Direct-to-Diffuse 
Signal Power Ratio = 0.01 
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B. Performance of Self-Normalization FFH/MFSK Receiver 

All parameters used in the graphs for the self-normalization receiver are the 
same as those used for the conventional receiver to make it easier to compare with 
each other as well as with the noise-normalization receiver. As we see from Fig. 4.17 
through Fig. 4.25, for a Rician fading channel, increasing diversity and increasing 
modulation order are both effective in minimizing the effects of worst case partial- 
band interference. 

Overall, the performance of the self-normalization receiver is almost the same 
as that of the noise-normalization receiver. There exists a small difference in the 
performance of the two receivers. This difference is expected due to the union bound 
used to obtain the FFH/MFSK self- normalized receiver performance. The union 
bound has the most pronounced effect on system performance when Eb/Nj is small 
and the modulation order is high. See, for example, Fig. 4.26 where the probability 
of bit error for a bit energy- to-jamming noise density ratio of 0 dB with a jamming 
ratio 1.0 is shown greater than unity. This is not valid for a probability. However, 
the error performance at low Eb/Nj is not important, and the error due to the union 
bound compared to the exact value at high Eb/Nj is not significant. 

The performance of the self-normalization receiver for a Rayleigh fading channel 
in worst case partial-band jamming is the same as for full band jamming (7 = 1) 
regardless of the parameters as shown in Fig. 4.27 through Fig. 4.29. Generally for 
a Rayleigh fading channel, diversity is effective in improving the performance of the 
self-normalization receiver. This is shown in Fig. 4.30. The effect of increasing the 
modulation order is less significant as shown in Fig. 4.26. As mentioned previously, 
the union bound makes the simulated performance worse than the true performance. 
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Figure 4.17: The Performance of the Self-Normalization FFH/MFSK non- 
coherent Receiver at Signal-to-Thermal Noise Ratio = 13.35 dB, Hopping 
Number per Symbol = 1, Order of Modulation = 4 and Direct-to-Diffuse 
Signal Power Ratio = 10 
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Self-Normalization FH/MFSK Receiver Performance 
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Figure 4.18: The Performance of the Self-Normalization FFH/MFSK non- 
coherent Receiver at Signal-to-Thermal Noise Ratio = 13.35 dB, Hopping 
Number per Symbol = 2, Order of Modulation = 4 and Direct-to-Diffuse 
Signal Power Ratio = 10 
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Figure 4 . 19 : The Performance of the Self-Normalization FFH/MFSK non- 
coherent Receiver at Signal-to-Thermal Noise Ratio = 13.35 dB, Hopping 
Number per Symbol = 3 , Order of Modulation = 4 and Direct-to-Diffuse 
Signal Power Ratio = 10 
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Self-Normal ization FH/MFSK Receiver Performance 
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Figure 4.20: The Performance of the Self-Normalization FFH/MFSK non- 
coherent Receiver at Signal-to-Thermal Noise Ratio = 13.35 dB, Hopping 
Number per Symbol = 4, Order of Modulation = 4 and Direct-to-Diffuse 
Signal Power Ratio = 10 
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Figure 4.21: The Performance of the Self-Normalization FFH/MFSK non- 
coherent Receiver at Signal-to-Thermal Noise Ratio = 13.35 dB, Hopping 
Number per Symbol = 5, Order of Modulation = 4 and Direct-to-Diffuse 
Signal Power Ratio = 10 
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Figure 4.22: The Performance of the Self-Normalization FFH/MFSK non- 
coherent Receiver at Signal-to-Thermal Noise Ratio = 13.35 dB, Hopping 
Number per Symbol = 6, Order of Modulation = 4 and Direct-to-Diffuse 
Signal Power Ratio = 10 
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Figure 4 . 23 : The Performance of the Self-Normalization FFH/MFSK non- 
coherent Receiver at Signal-to-Thermal Noise Ratio = 13.35 dB, Hopping 
Number per Symbol = 4 , Order of Modulation = 8 and Direct-to-Diffuse 
Signal Power Ratio = 10 
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Figure 4.24: The Performance of the Self-Normalization FFH/MFSK non- 
coherent Receiver at Signal-to-Thermal Noise Ratio = 13.35 dB, Hopping 
Number per Symbol = 4, Order of Modulation = 16 and Direct-to-Diffuse 
Signal Power Ratio = 10 
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Figure 4.25: The Performance of the Self-Normalization FFH/MFSK non- 
coherent Receiver at Signal-to-Thermal Noise Ratio = 13.35 dB, Hopping 
Number per Symbol = 4, Order of Modulation = 32 and Direct-to-Diffuse 
Signal Power Ratio = 10 
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Figure 4.26: The Performance of the Self-Normalization FFH/MFSK non- 
coherent Receiver at Signal-to-Thermal Noise Ratio = 20.0dB, Hopping 
Number per Symbol = 4, Jamming Ratio = 1 and Direct-to-Diffuse Signal 
Power Ratio = 0.01 
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Figure 4.27: The Performance of the Self-Normalization FFH/MFSK non- 
coherent Receiver at Signal-to-Thermal Noise Ratio = 20.0 dB, Hopping 
Number = 1, Order of Modulation = 4 and Direct-to-Diffuse Signal Power 
Ratio = 0.01 
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Figure 4.28: The Performance of the Self-Normalization FFH/MFSK non- 
coherent Receiver at Signal-to-Thermal Noise Ratio = 20.0dB, Hopping 

Number = 4, Order of Modulation = 4 and Direct-to-Diffuse Signal Power 
Ratio = 0.01 



58 



o 




CO 

o 

>- 



JD 

no 

_□ 

O 

L_ 

CL_ 



Self-Normalization FH/MFSK Receiver Performance 
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Figure 4.29: The Performance of the Self-Normalization FFH/MFSK non- 
coherent Receiver at Signal-to-Thermal Noise Ratio = 20.0dB, Hopping 
Number = 6, Order of Modulation = 32 and Direct-to-Diffuse Signal Power 
Ratio = 0.01 
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Figure 4.30: The Performance of the Self-Normalization FFH/MFSK non- 
coherent Receiver at Signal-to-Thermal Noise Ratio = 20.0dB, Jamming 

Ratio 1, Order of Modulation = 8 and Direct-to-Diffuse Signal Power 
Ratio = 0.01 
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V. CONCLUSION 



The probability of bit error performance for a conventional noncoherent FFH/MFSK 
receiver has been obtained for Rayleigh- and Rician-faded channel with partial- 
band jamming. The effect of fading is detrimental for both fading models, but fast 
frequency-hopping in general provides a means to overcome, at least partially, fading 
effects. However, the conventional receiver is severely affected by worst case partial- 
band jamming. It does not show any improvement in worst case performance with 
high diversity or high order of modulation. The probability of bit error of the conven- 
tional receiver in worst case partial-band jamming is always between 10~ 2 and 10 -3 
at Eb/Nj of 20 dB regardless of E^/No, the order of modulation, diversity factor, or 
direct-to diffuse power ratio. The probability of bit error of the self-normalization re- 
ceiver in worst case partial-band jamming is the same as that of conventional receiver 
for slow frequency-hopping (L=l). However, for fast frequency-hopping, the perfor- 
mance becomes better and better with increasing diversity factors: the probability of 
bit error is 5 x 10“ 4 at L=2, 2 x 10“ 4 at L=3, 10“ 4 at L=4 and 8 x 10“ 5 at L=5 and 
L=6. When the diversity factor is 6, the worst case performance is almost same as 
that of full band jamming case (jamming ratio 1.0). At this point the jammer does 
not need any information about the communicator to optimize the jamming effect by 
managing the jamming power. 

The self-normalization receiver with diversity provides very good immunity to 
worst case partial-band jamming and has a performance similar to that of the noise- 
normalization receiver. The order of modulation also improves the overall perfor- 
mance of the self-normalization receiver, but partial-band jamming affects higher 
orders of modulation more than the lower ones. Thus, if the communicator uses a 
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higher order of modulation and wants to suppress all partial-band jamming effects, 
he should use a higher hopping rate. 
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APPENDIX A 



% This is the Matlab program to produce the probability of bit % 
error of 

% Conventional FH/MFSK quadratic receiver. 

% Lee, Kang Yeun (Major ROKAF) 

clear 

% 3 signal-to-noise ratios to be implemented (13.35 db, 16 db and 
% 20 db) 

SNR= [21. 62 7 18 52 4 DO 39 . 81071706DO 100. 0D0] ; 

% 2 direct-to-dif fuse signal power ratios (0.01 and 10) 

DDR= [ 0 . 01 DO 10.0D0]; 

% 4 jamming ratio (1, 0.1, 0.01 and 0.001) 

JRATIO= [ 1 . 0D0 0.1D0 0.01D0 0.001D0]; 
data=zeros (8 , 51) ; 

% set the thermal plus jamming noise power 1 
rhokl=l ; 

% basic sequence length is 1024 
lim=1024 ; 
n=2*lim; 

zpad=zeros (1 , lim) ; 
opad=ones (1 , lim) ; 

% hopping number (1 to 6) 
for L=l:6, 

for jl=l:2, psi=DDR(jl) ; 
for j2=l:3, nu=SNR(j2); 
for j3=l:4, gama=JRATIO ( j 3) ; 

% the order of modulation (M-ary : 4, 8, 16, 32-ary) 
for j4=2 : 5 ; M=2 A j4; 

for j5=0: 2 : 50, iota=10. 0D0 A (0. lD0*j5) ; 

% diffuse signal power ( rho ) 

rho=nu*gama*iota/ (nu+gama*iota) / (psi+1) * j 4/L; 

% direct signal power ( alpha ) 
alpha=rho*psi ; 

% thermal noise power ( rhok2 ) 
rhok2=iota*gama/ ( iota*gama+nu) ; 

% set the sequence limit roughly based on the singal-to-noise ratio 
% (nu) 

limx=80*nu ; 

% divide the range by 100 to check the proper limit 
delta=limx/100 ; 

% select almost the middle value of the increment to compute the 
% sequence 

% of the probability density function. 
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x=0 . 49*delta : delta : limx-0 . 51*delta ; 

g=exp (- (x+2*L*alpha) /2/ (rhokl+rho) ) /2/ (rhokl+rho) . ... 

* x/2/L/alpha) . A . . . 

( (L-l) /2) . *abs ( j A (L-l) *bessel (L-l, j*sqrt (2*L*alpha*x) . . . 

/ (rhokl+rho) ) ) *delta; 

% check the sequence whether it has invalid sequence or not 
fi=find(isnan(g) ) ; 
g ( f i ) =zeros ( 1 , length ( f i ) ) ; 

% if the sequence has only Os, the limit is too wide. Thus set the 
% first 

% increment as the new limit 
if all (g==0)==l, 
gmin=l ; 

% if the sequence is valid, check the point around 1.0e-6 to select 
% the 

% new limit 

else, [xx,gmax]=max(g) ; 
g(l:gmax)=ones(l,gmax) *xx; 

[xx, gmin]=min ( (g-1.0e-6) . A 2) ; 
end 

1 imx=gmin*delta ; 

% if the new limit is first or second value of the increment 
% redo above procedure 
while gmin <= 2 , 

delta=l imx/100 ; x=0 . 5*delta : delta : limx-0 . 5*delta ; 
g=exp(-(x+2*L*alpha)/2/( rhokl+rho) )/2/ . . . 

(rhokl+rho) .* (x/2/L/alpha) . A ... 

( (L-l) /2) . *abs ( j A (L-l) *bessel (L-l , j *sqrt (2*L*alpha*x) . . . 

/ (rhokl+rho) ) ) *delta ; 
fi=f ind ( isnan (g) ) ; 
g (fi)=zeros(l, length (fi) ) ; 
if all (g==0) ==1, 
gmin=l ; 

else, [xx,gmax]=max(g) ; 
g(l:gmax)=ones(l,gmax) *xx; 

[xx, gmin]=min ( (g-1.0e-6) . A 2) ; 
end 

1 imx=gmin*delta ; 
end 

% after setting the limit, whole range is divided by 1024 
% to produce the sequence 

delta=limx/lim;x=0. 49 *delta: delta: limx-0. 51 *delta; 

% if jamming ratio is 1, the computation is simpler than others 
if j 3 == 1, 

gs=- (x+2*L*alpha)/2/ (rhokl+rho) -log (2* (rhokl+rho) ) 
+log(x/2/L/alpha) * ... 

(L-l) /2+log (abs ( j A (L-l) *bessel (L-l, j *sqrt (2*L*alpha*x) . . . 

/ (rhokl+rho) ) ) ) ; 
gs=exp(gs) *delta; 
f i=f ind ( isnan (gs) ) ; 
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gs (f i) =zeros (1, length ( f i) ) ; 
if all (gs==0)==l, 
gs=zeros ( 1, lim) ; 
gs(l)=l; 

else,gs=gs/sum(gs) ; 
end 

ps=log ( 1/2/rhokl) *L+log(x) * (L-l) -log(fact(L-l) ) - x/2/rhokl; 

ps=exp (ps) *delta ; 

f i=f ind ( isnan (ps) ) ; 

ps (f i) =zeros (1, length ( f i) ) ; 

if all (ps==0) ==1 , 

ps=ones ( 1 , lim) ; 

else,ps=ps/sum(ps) ; 

ps=cumsum (ps) ; 

ps=ps/max (ps) ; 

end 

data (count, j 5+1) =M/2/ (M-l) * ( l-sum(gs . *ps . A (M-l) ) ) ; 
else 

% if jamming ratio is not 1 , the hopping number has to be 
% considered 
for i = 0 : L 
% if no hop is jammed 
if i == 0, 
delta=limx/100; 

x=0 . 5*delta : delta : limx-0 . 5*delta ; 

g=- (x+2*L*alpha) /2/ (rhok2+rho) -log(2* (rhok2+rho) ) ... 

+log (x/2/L/alpha) * ... 

(L-l) /2+log(abs ( j A (L-l) *bessel (L-l, j*sqrt(2*L*alpha*x) . . . 

/ (rhok2+rho) ) ) ) ? 

g=exp (g) *delta ; 

fi=find(isnan(g) ) ; 

g ( f i) =zeros ( 1 , length ( f i) ) ; 

if all (g==0) ==1 , 

gmin=l ; 

else, [xx, gmax] =max (g) ;g(l:gmax)=ones(l,gmax) *xx; 

[xx, gmin]=min ( (g-1. 0e-6) . A 2) ; 
end 

vlimx=gmin*delta ; 

% check the sequence limit again 
while gmin <= 2, 
delta=vlimx/100 ; 

x=0 . 5*delta : delta : vlimx-0 . 5*delta ; 

g=- (x+2*L*alpha)/2/ (rhok2+rho) -log(2* (rhok2+rho) ) ... 

+log(x/2/L/alpha) * ... 

(L-l) /2+log (abs ( j A (L-l) *bessel (L-l , j *sqrt (2*L*alpha*x) . . . 

/ ( rhok2+rho) ) ) ) ; 
g=exp (g) *delta ; 
fi=find(isnan(g) ) ; 
g ( f i) =zeros ( 1 , length ( f i) ) ; 
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if all (g==0) ==l, 

gmin=l ; 

else, 

[xx,gmax]=max(g) ;g(l:gmax)=ones(l,ginax) *xx; 

[xx,gmin]=min( (g-1.0e-6) . A 2) ; 
end 

vl imx=gmin*del ta ; 
end 

delta=vlimx/lim; 

x=0 . 49*delta : delta : vlimx-0 . 51*delta ; 

gs=- (x+2*L*alpha) /2/ (rhok2+rho) -log(2* (rhok2+rho) ) ... 

+log(x/2/L/alpha) * ... 

(L-l)/2+log(abs(j A (L-l) *bessel (L-l, j*sgrt(2*L*alpha*x) . . . 

/ (rhok2+rho) ) ) ) ; 

gs=exp (gs) *delta ; 

fi=find(isnan(gs) ) ; 

gs (fi)=zeros(l, length (fi) ) ; 

if all (gs==0) ==1, 

gs=zeros (1, lim) ; 

gs(l)=l; 

else, 

gs=gs/sum(gs) ; 
end 

ps=-log(2*rhok2) *L+log(x) * (L-l) -log(fact(L-l) ) - x/2/rhok2; 

ps=exp (ps) *delta ; 

fi=find(isnan(ps) ) ; 

ps (f i)=zeros(l, length (fi) ) ; 

if all (ps==0) ==1, 

ps=ones(l,lim) ; 

else, ps=ps/sum(ps) ; 

ps=ps/sum(ps) ; 

ps=cumsum ( ps ) ; 

ps=ps/max (ps) ; 

end 

% if all hops are jammed 
elseif i == L, 
delta=limx/100 ; 

x=0. 5*delta: delta: limx-O. 5*delta; 

g=- (x+2*i*alpha) /2/ (rhokl+rho) -log(2* (rhokl+rho) ) ... 

+log(x/2/ i/alpha) * ... 

(i-l)/2+log(abs( j A (i-1) *bessel (i-1, j*sqrt(2*i*alpha*x) . . . 

/ (rhokl+rho) ) ) ) ; 

g=exp (g) *delta ; 

f i=f ind ( isnan (g) ) ; 

g(fi)=zeros(l, length (fi) ) ; 

if all (g==0)==l, 

gmin=l ; 

else, 

[xx, gmax]=max (g) ;g(l:gmax)=ones(l,gmax) *xx? 

[xx, gmin]=min ( (g-1.0e-6) . A 2) ; 
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end 

ulimx=gmin*delta ; 

% check the sequence limit again 
while gmin <=2, 
delta=ulimx/100 ; 

x=0 . 5*delta : delta : ul imx-0 . 5*delta ; 

g=- (x+2*i*alpha) /2/ (rhokl+rho) -log (2* (rhokl+rho) ) ... 

+log (x/2/i/alpha) * ... 

(i-l)/2+log(abs( j A (i-1) *bessel ( i— 1 , j*sqrt (2*i*alpha*x) . . . 

/ (rhokl+rho) ) ) ) ; 
g=exp (g) *delta ; 
f i=f ind (isnan (g) ) ; 
g ( f i ) =zeros ( 1 , length ( f i ) ) ; 
if all (g==0)==l, 
gmin=l ; 
else, 

[xx, gmax]=max(g) ; 
g(l:gmax)=ones(l,gmax) *xx; 

[xx, gmin]=min ( (g-1 . Oe-6) . A 2 ) ; 
end 

ulimx=gmin*delta ; 
end 

delta=ulimx/lim; 

x=0 . 49*delta : delta : ul imx-0 . 51*delta ; 

gs=-(x+2*i*alpha)/2/ (rhokl+rho) -log (2* (rhokl+rho) ) ... 

+log (x/2/i/alpha) * ... 

( i— 1 ) /2+log(abs ( j A ( i — 1 ) *bessel (i-1, j *sqrt (2*i*alpha*x) . . . 

/ (rhokl+rho) ) ) ) ; 
gs=exp(gs) *delta; 
f i=find (isnan(gs) ) ; 
gs (fi)=zeros(l, length (fi) ) ; 
if all (gs==0) ==1, 
gs=zeros (1 , lim) ; 
gs ( 1) =1 ; 
else, 

gs=gs/sum(gs) ; 
end 

ps=-log (2*rhokl) *i+log (x) *( i-1) -log ( fact (i-1) ) - x/2/rhokl; 

ps=exp(ps) *delta; 

if all (ps==0) ==1, 

ps= ones ( l , lim) ; 

else, 

fi= f ind ( isnan (ps) ) ; 

ps ( f i) =zeros ( l, length ( f i ) ) ; 

ps = cumsum(ps) ; 

ps = ps/max (ps); 

end 
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else 

% if i hops are jammed and L-i hops are not jammed 
% check the sequence limit again 
rhok=rhokl ; 
delta=limx/100? 

x=0 . 5*delta : delta : limx-0 . 5*delta ; 

g=-(x+2*i*alpha)/2/(rhok+rho)-log(2*(rhok+rho) ) ... 

+log (x/2/i/alpha) * ... 

(i-1) /2+log(abs ( j A (i-1) *bessel (i-1, j*sqrt (2*i*alpha*x) 

/ (rhok+rho) ) ) ) ; 

g=exp (g) *delta ; 

fi=find(isnan(g) ) ; 

g(fi)=zeros (1, length (fi) ) ; 

if all (g==0) ==1, 

gmin=l ; 

else, 

[xx,gmax]=max(g) ; 
g ( 1 : gmax) =ones ( 1 , gmax) *xx ; 

[xx,gmin]=min( (g-1.0e-6) . A 2) ; 
end 



wlimx=gmin*delta; 
while gmin <=2, 
delta=wlimx/100 ; 

x=0. 5*delta:delta:wlimx-0. 5*delta; 

g=- (x+2*i*alpha) /2/ (rhok+rho) -log (2* (rhok+rho) ) ... 
+log (x/2/i/alpha) * ... 

(i-l)/2+log(abs( j A (i-1) *bessel (i-1, j*sqrt (2*i*alpha*x) 

/ (rhok+rho) ) ) ) ? 

g=exp (g) *delta ; 

fi=find(isnan(g) ) ; 

g(fi)=zeros(l,length(fi) ) ; 

if all (g==0) ==1, 

gmin=l ; 

else, 

[xx,gmax]=max(g) ; 
g ( 1 : gmax) =ones ( 1 , gmax) *xx ; 

[xx,gmin]=min( (g-1.0e-6) . A 2) ; 
end 

wl imx=gmin*delta ; 
end 

rhok=rhok2 ; 
delta=limx/100 ; 

x=0. 5*delta:delta: limx-O. 5*delta; 

g=- (x+2* (L-i) *alpha)/2/ (rhok+rho) -log (2* (rhok+rho) ) 
+log(x/2/ (L-i) /alpha) * ... 

(L-i-l)/2+log(abs ( j A (L-i-1) *bessel (L-i-1, j *sqrt (2* (L-i) 
*alpha*x)/ ... 

(rhok+rho) ) ) ) ? 
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g=exp (g) *delta ; 

f i=f ind ( isnan (g) ) ; 

g(fi)=zeros(l,length(fi) ) ; 

if all (g==0) ==1, 

gmin=l ; 

else, 

[xx,gmax]=max(g) ; 
g(l:gmax)=ones(l,gmax) *xx; 

[xx,gmin]=min( (g-1. Oe-6) . A 2) ; 
end 

zlimx=gmin*delta; 

% check the sequence limit again 
while gmin <=2, 
delta=zlimx/100 ; 

x=0. 5*delta: delta: zlimx-O. 5*delta; 

g=- (x+2* (L-i) *alpha) /2/ (rhok+rho) -log (2* (rhok+rho) ) ... 
+log (x/2/ (L-i) /alpha) * ... 

(L-i-1) /2+log (abs (j A (L-i-1) *bessel (L-i-1 , j*sqrt (2* (L-i) 

*alpha*x) / (rhok+rho) ) ) ) ; 

g=exp (g) *delta ; 

fi=f ind ( isnan (g) ) ; 

g(fi)=zeros (1, length (fi) ) ; 

if all (g==0) ==1 , 

gmin=l ; 

else, 

[xx,gmax]=max(g) ; 
g(l:gmax)=ones (l,gmax) *xx; 

[xx, gmin]=min ( (g-1. Oe-6) . A 2) ; 
end 

zlimx=gmin*delta; 

end 

% produce the sequency of jammed hop of signal branch 
del ta=wl imx/1 im ; 

x=0. 49 *delta: delta :wlimx-0. 51*delta ; 

gl=- (x+2*i*alpha) /2/ (rhokl+rho) -log (2* (rhokl+rho) ) ... 

+log (x/2/i/alpha) * ... 

( i-1 ) /2+log (abs(j A (i-1) *bessel (i-1, j*sqrt(2*i*alpha*x) . 
/ (rhokl+rho) ) ) ) ; 
gi =exp (gl) *delta ; 
f i=f ind (isnan(gl) ) ; 
gl (f i) =zeros (1, length ( f i) ) ; 
if all (gl==0) ==1, 
gl=zeros (1, lim) ; 
gl(l)=l; 
else, 

gl=g 1/sum (gl) ; 
end 

% produce the sequence of nonjammed hop of signal branch 
if wlimx >= 100*zlimx 
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gs=[gl zpad] ; 
else 

g2=- (x+2* (L-i) *alpha)/2/ (rhok2+rho) -log(2* (rhok2+rho) ) ... 

+log (x/2/ (L-i) ... 

/alpha) * (L-i-1) /2+log(abs (j A (L-i-1) *bessel (L-i-1 , ... 

j*sqrt (2* (L-i) ... 

*alpha*x) / (rhok2+rho) ) ) ) ; 
g2=exp(g2) *delta; 
f i=f ind (isnan (g2) ) ; 
g2 (fi)=zeros(l, length (fi) ) ? 
if all (g2==0)==l, 
g2=zeros (1, lim) ; 
g2 (l)=l; 
else, 

g2=g2/sum(g2) ; 
end 

% produce the sequence of signal branch 
gs = abs(ifft(fft(gl,n) .*fft(g2,n) ,n) ) ; 
gs = gs/sum(gs) ; 
end 

% produce the sequence of jammed hop of noise branch 

pl=-log (2*rhokl) *i+log(x) * (i-1) -log (fact (i-1) ) -x/2/rhokl ; 

pl=exp (pi) *delta ; 

fi=find(isnan(pl) ) ; 

pi (fi) =zeros (1, length (fi) ) ; 

if all (pl==0)==l, 

pl=zeros (1, lim) ; 

pl(l)=l; 

else, 

pl=pl/sum(pl) ; 
end 

% produce the sequence of non jammed hop of noise branch 

p2=-log(2*rhok2) * (L-i)+log(x) * (L-i-1) -log(fact (L-i-1) ) -x/2/rhok2 

p2=exp (p2 ) *delta ; 

fi=find(isnan(p2) ) ; 

p2 (fi)=zeros(l, length (fi) ) ; 

if all (p2==0) ==1 , 

p2=zeros (1 , lim) ; 

p2 ( 1 ) =1 ; 

else, 

p2=p2/sum(p2) ; 
end 

% produce the sequence of noise branch 
ps = abs(ifft (fft (pl,n) .*fft(p2,n) ,n) ) ; 
ps = cumsum(ps) ; 
ps = ps/max (ps) ; 
end 
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% compute the probability of bit error using sequence 
tdata ( j 5+1) =tdata ( j 5+1) +bino (L, i) *gama A i* ( 1-gama) A (L-i) 
*M/2* (l-sum(gs. *ps. A (M-l) ) ) ; 
end, end, end 

data=[data;td] ; 

save datafile data; 

end , end , end , end , end 
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APPENDIX B 



% This is the program to produce the probability of bit error of 
% Self-Normalization FH/MFSK Receiver 
% Lee, Kang Yeun (Major, ROKAF) 

% three special routines 

% fact : factorial routine 

% bino : binomial coefficient routine 

% ereval : error probability computation routine using fft 

% 3 signal-to-noise ratio(13 . 35db, 16db, 20db) 

SNR=[21 . 62 7 18 52 4 DO 39 . 81071706D0 100. 0D0] ; 

% 3 direct-to-dif fuse signal power ratio (0.01, 1, 10) 

DDR= [ 0 . 01D0 1.0D0 10.0D0]; 

% 4 jamming ratio 

JRATI01= [ 1 . 0D0 0.2D0 0.1D0 0.05D0 0.01D0]; 

JRATIO2=[1.0D0 0.2D0 0.05D0 0.01D0 0.001D0]; 

% Bandwidth-to-Symbol Rate Ratio (BW/Rs : 1 2 ) 

mu=l . 0D0 ; 

% M-ary 
M=16 . 0D0 ; 

% Multichannel factor w=logM/log2 
omega=log(M) /log (2 . 0D0) ; 

% Bit SNR 13.35 db 
nu=SNR(l) ; 

count=l ; 

% Hopping number (1 to 6) 
for Jl=l : 6 
L=J1 ; 

% Direct-to-Dif fuse signal ratio (0.01 1 10 ) 

for J3=l:3 
psi=DDR(J3) ; 

% Direct SNR of NOJAMMING CASE 

rho=omega*nu*psi/mu/ (psi+1. 0D0)/L; 

% Diffuse SNR of NOJAMMING CASE 
xi=omega*nu/mu/ (psi+1 . 0D0) /L; 

% produce the sequence of nojamming case 
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z=-0.999D0:0.002D0:0.999D0; 

NOJAM=exp (rho* ( z-1 . 0D0) ./ (2 . ODO+xi* ( 1 . ODO-z ) ) ) . ... 

* (1. ODO+z) . A (M-2 . ODO) . . . 

./ (2. ODO+xi* (1. ODO-Z) ) . *( ( 2. ODO* (1. ODO+xi) ) ./ ... 
(xi* ( 1. ODO-z)+2 . ODO) . . . 

. * (l.ODO+rho* (1. ODO+z) ./( (2. ODO+xi* (1. ODO-z) ) ... 

*(1. ODO+xi) ) ) +M-2 . ODO) ; 



z=0.001D0:0.002D0:0.999D0; 

B=rho*z ./ ( (1.0D0+xi*(l. ODO-z) )*(l. ODO+xi) ) ; 

E=exp(rho* (z-1. ODO) ./ (1. ODO+xi* (1. ODO-z) ) ) ; 

TT=zeros (M-l , 500) ; 
for N=0 : M-2 

TT (N+l , : ) =fact (M-2 . ODO) /fact (M-2 . ODO-N) /fact (N) A 2*B. A N 7 
end 

Tl=sum(TT) .*( (1. ODO-z) *(1. ODO+xi) ) . A (M-2) .*z. ... 

* (M-2. ODO) ./(l. ODO+z) . . . 

. /(1.0D0+xi*(l. ODO-z) ) . A (M-l ) . *E ; 

TT=zeros (M-2 , 500) ? 
for N=0 :M-3 
TTl=zeros (M-l-N, 500) ; 
for IQ=0 : M-2-N 

TT1 (IQ+1, : ) =fact (M-2 . ODO-N) /fact (M-2 . ODO-N-IQ)/fact (IQ) A 2*B. A IQ 
end 

TEMP=2 . ODO* ( (1. ODO-z) * (1. ODO+xi) ) . A (M-3-N) . * (1. ODO+xi) . *z. A 2 
.*(1. ODO+Z) . A (N-1) . * (M-2. ODO-N) 7 
TEMP=TEMP. /(I. ODO+xi* (1. ODO-Z) ) . A (M-l . ODO-N) . *E . *sum (TT1) ; 
TTl=zeros (M-2-N ,500) ? 
for IQ=0 :M-3-N 

TT1 ( IQ+1 , : )=f act (M-3. ODO-N) /fact (M-3 . ODO-N-IQ) /fact (IQ) A 2*B. A IQ 
end 

TEMP1=( (1. ODO-z) *(1. ODO+xi) ) . A (M-3-N) . *z .*( 1 . ODO+Z ) . . . 

. A N.* (M-2. ODO) ./(1.0D0+xi*(l. ODO-Z) ) . A (M-2-N) .*E; 
if IQ == 0, 

TT (N+l , : )=TEMP+TEMP1.*TT1; 
else 

TT (N+l , : )=TEMP+TEMPl.*sum(TTl) ? 

end 

end 

Tl=Tl+sum (TT) ; 

% sequence of nojanuning case 
NOJAM (501: 1000) =NOJAM( 501: 1000) -Tl; 

Ps=zeros(5, 41) ; 



% jamining ratio 
for J4=l: 5 
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if J3==3 

gamma=JRATI02 ( J4 ) ; 
else 

gamma=JRATI01 (J4 ) ; 
end 

% Signal-to-Jamming noise ratio ( 0 db - 40 db ) 
for J5=0:40 

iota=10. ODO A (0. 1D0*J5) ; 

% Direct SNR of JAMMING case 

rho=omega*nu*iota*psi*gamma/mu/ (psi+1 . 0D0)/ (iota*gamina+nu) /L; 

% Diffuse SNR of JAMMING case 

xi=omega*nu*iota*gamma/mu/ (psi+1. 0D0) / (iota*gamma+nu)/L; 

% produce the sequence when the signal is jammed by noise jamming 

z=-0.999D0:0.002D0:0.999D0; 

JAM=exp(rho*(z-1.0D0) ./(2.0D0+xi*(l. ODO-z) ) ) . ... 

* (1 . ODO+z ) . A (M-2 . ODO) . . . 

./ (2 . ODO+xi* (1 . ODO-z) ) . ... 

*( (2.0D0*(1. ODO+xi) ) ./(xi*(l. ODO-Z) +2 . 0D0) . . . 

. *(1.0D0+rho*(l. ODO+z) ./( (2. ODO+xi* (1. ODO-z ) ) . . . 

*(1. ODO+xi) ) ) +M-2 . ODO) ; 



Z=0 . 001D0 : 0 . 002 DO : 0 . 999D0 ; 

B=rho*z ./ ( (1. ODO+xi* (1 . ODO-z) ) * (1 . ODO+xi) ) ; 

E=exp ( rho* ( z-1 . ODO ) ./ (1 . ODO+xi* (1 . ODO-z) ) ) ; 

TT=zeros(M-l, 500) ; 
for N=0 :M-2 

TT (N+l , : )=fact (M-2. ODO) /fact (M-2. ODO-N) /fact (N) A 2*B. A N; 
end 

Tl=sum(TT) .*( (1. ODO-Z) * (1. ODO+xi) ) . A (M-2) . ... 

*z.* (M-2. ODO) ./(l. ODO+z) . . . 

./(I. ODO+xi* (1. ODO-Z) ) . A (M-l) . *E ; 

TT=zeros (M-2 , 500) ; 
for N=0 : M-3 
TTl=zeros (M-l-N, 500) ; 
for IQ=0 :M-2-N 

TT1 (IQ+1, : ) =fact (M-2 . ODO-N) /fact (M-2 . ODO-N-IQ) /fact (IQ) A 2*B. A IQ; 
end 

TEMP=2. ODO* ((1. ODO-Z) *(1. ODO+xi)) . A (M-3-N) .*(1. ODO+xi) .*Z. A 2 ... 

.*(1. ODO+z) . A (N-l ) .* (M-2. ODO-N) ; 

TEMP=TEMP. / ( 1 . ODO+xi* ( 1 . ODO-z ) ) . A (M-l . ODO-N) . *E . *sum (TT1) ; 
TTl=zeros (M-2-N, 500) ; 
for IQ=0 : M-3 -N 

TT1 (IQ+1 , : ) =fact (M-3 . ODO-N) /fact (M-3 . ODO-N-IQ) /fact (IQ) A 2*B. A IQ; 
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end 

TEMP1= ( ( 1 . ODO-Z ) * ( 1 . ODO+xi) ) . A (M-3-N) . *z . * (1 . ODO+z) . . . 

. A N. * (M-2 . ODO) ./(!• ODO+xi* (1. ODO-Z) ) . A (M-2-N) . *E; 
if IQ == 0, 

TT (N+l , : ) =TEMP+TEMP1 . *TT1 ; 
else 

TT (N+l , : ) =TEMP+TEMP1 . * sum ( TT1 ) ; 

end 

end 

Tl=Tl+sum(TT) ; 

% the sequence of jamming case 
JAM (501: 1000) =JAM( 501: 1000) -Tl; 



% 'bino' and 'ereval' are the functions for binomial and error 
% evaluation. 



% Ps ( i) <= (M-l) * Pr ( z < 0 ) 
% Pb = (M/2)/ (M-l) *Ps 
% Pb <= (M/2) * Pr( z < 0 ) 



% probability of symbol error 
for i=0:L 

Ps (J4 , J5+1) =Ps (J4 , J5+1) +bino (L, i) *gamma A i* (1-gamma) A (L-i) . . . 

*ereval (NOJAM, L-i, JAM, i) ; 



end 



% probability of bit error 

Ps (J4 , J5+1) = (M/2 . ODO) *Ps ( J 4 , J5+1) ; 

end 

end 



% Save the Results and Plot 



save data Ps 



end 

end 
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